Transport, optical and electronic properties of the half metal CrC>2 
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The electronic structure of CrC>2 is critically discussed in terms of the relation of existing ex- 
perimental data and well converged LSDA and GGA calculations of the electronic structure and 
transport properties of this half metal magnet, with a particular emphasis on optical properties. We 
find only moderate manifestations of many body effects. Renormalization of the density of states 
is not large and is in the typical for transition metals range. We find substantial deviations from 
Drude behavior in the far-infrared optical conductivity. These appear because of the unusually low 
energy of interband optical transitions. The calculated mass renormalization is found to be rather 
sensitive to the exchange-correlation functional used and varies from 10% (LSDA) to 90% (GGA), 
using the latest specific heat data. We also find that dressing of the electrons by spin fluctuations, 
because of their high energy, renormalizes the interband optical transition at as high as 4 eV by 
about 20%. Although we find no clear indications of strong correlations of the Hubbard type, strong 
electron-magnon scattering related to the half metallic band structure is present and this leads to 
a nontrivial temperature dependence of the resistivity and some renormalization of the electron 
spectra. 



C1O2 has been-, attracting substantial theoreticalErB 
and experimentaluilil interest in recent years, in part 
because of its practical importance in magnetic record- 
ing, and in part due to fundamental interest in its half- 
metallic electronic structure and the implications of this. 
In half metals, one spin channel is metallic and the other 
is semiconducting. This results in unusual transport 
properties, including anomalously Large, but still metal- 
lic, resistivities at high temperatureEj, coupled with high 
mobilities at low temperature. Another consequence of 
the half metal behavior is nontrudaL magnetoresistance 
due to spin-dependent tunnelingEBiij, with various po- 
tential applications. 

Despite a large number of published band structure 
calculations, no full potential, all-electron calculations 
have been reported fot_Cr02 except for calculations of 
Sorantin and SchwarzE-3, done in a study in chemical 
bonding of the rutile structure oxides. The calculations 
of Refs. used atomic sphere approximations for 

the crystal potential, while those of Ref. ^ used a muffin- 
tin approximation. Note that nonspherical part of the 
potential is of particular significance for a material with 
low site symmetries like rutile Cr02- The only nonsphc 
ical calculations reported so far are those of Lewis et ai 
however the accuracy of these calculations was limited 
by the use of a pseudopotential based method. 

As was observed early oni, the exchange splitting is 
very sensitive to computational details, hence the ne- 
cessity of well-converged, all-electron, and full-potential, 
state-of-the-art band structure calculations to establish 



the predictions of the LSDA which may then be compared 
with experiment to clarify the possible role of many body 
effects beyond the mean field. Currently even such a ba- 
sic quantity as density of states at the Fermi level, N(0), 
varies among different calculations from 0.69 states/Ry 
f.u. (Ref. g) to 2.34 states/f.u.(Ref. |). Given that the 
reported experimental low temperature specific heat co- 
efficients range from 2.5 mJ/K 2 (Ref. g) to 7 mJ/K 2 (Ref. 
|l5| ), such a situation creates confusion about the physi- 
cal nature of this compound, including speculations that 
strong correlations may lead to unusual renormalizations 
in Cr0 2 . 

As mentioned, here we report well converged full- 
potential all-electron band structure calculations along 
with an analysis of the latest experimental data. Our 
conclusion is that the existing experimental data are rec- 
oncilable with the calculated LSDA electronic structure 
without additional strong correlation effects, but that 
strong electron-magnon scattering (as expected in half 
metals) is presents and renormalizes the electron bands. 
With such renormalization the LSDA electronic struc- 
ture suffice to understand reasonably well the body of 
experimental results. 

The calculations were performed using the gen- 
eral potential linearized augmented plane wavp_mcthod 
(LAPW) , which is described in detail elsewheretSEZI. Lo- 
cal orbital extensions were used to relax linearization er- 
rors and to treat semicore states. Well converged basis 
sets of over 600 functions for the six atom unit cell were 
used. The radii of the muffin-tin spheres were 2 a.u. for 
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Cr and 1.55 a.u. for O. The Hedin-Lundquist exchange- 
correlation potential was used for most calculations (but 
see below). The Brillouin zone averages during the it- 
erations to self consistency were performed with 105 in- 
equivalent k-points, while the zone integrals for calculat- 
ing transport and related properties utilized the tetra- 
hedron method with at least 770 inequivalent k-points. 
Convergence was carefully checked. 

The resulting band structure and electronic density of 
states (DOS) are shown in Fig.[l] and ^. The DQS,agrees 
well with that reported by Sorantin and SchwarzEJ. First 
of all, let us note that the LSDA bandwidth W is of 
the order of 10 eV. This is to be compared with the 
Hubbard repulsion U ~ 3 eV from constrained LSDA 
calculations!! The ratio of U /W w 0.3 is much smaller 
than in strongly correlated oxides of Ni, Fe or Cu. One 
can expect at best a moderate effect of electronpelectron 
correlations beyond the mean-field LSDA result£3. 

In this connection, it is instructive toJook at the recent 
LSDA+U calculations of Korotin et cm. They compared 
their results with an earlier LSDA band structured, and 
found some differences, mostly amounting to an upward 
shift of the unoccupied ci-states by « 2 eV. However the 
reported LSDA calculations deviate among themselves 
by a comparable amount, and it worth trying to sort 
out the differences genuinely related with the Hubbard 
U from those related to the computational details. 

The results of Korotin et afl mainly differ from those 
of Schwarzl in the position of the unoccupied spin-down 
Ci-d band. Compared to the present LSDA results (Figs. 
|l|,||) this band is pushed up by about 1.4 eV (raflre than 
twice less than U, due to the Cr-0 hybridizationlij) . Thus 
the Fermi level appears to be somewhat lower than in the 
middle of the gap, while in the LAPW calculations re- 
ported here, as well as in Ref. [I], it is closer to the top of 
the gap. The question is to what extent this difference 
reflects real physics, that is, renormalization of the bands 
due to electron-electron correlations beyond LSDA, and 
to what extent it reflects the difference in the computa- 
tional technique and associated approximations. If the 
latter dominates, full potential calculations are preferred 
over Atomic Sphere Approximation (ASA) calculations 
as in Ref. 0. On the other hand, if the correlations beyond 
the LSDA are very strong, the LDA+U calculations may 
be superior regardless of the spherical approximation. It 
isrworth noting that although the statement of Korotin et 
aH, that LSDA+U "is superior to LSDA since it can in- 
deed yield a gap if the local Coulomb interaction is large 
enough", is often true, it is not so clear when the local 
interaction is small compared to the band width, and the 
material is actually a metal. In fact, when U/W ~ 0.3, 
as in Cr02, one expects the crude way in which LSDA+U 
includes correlations to be not better (if not worse) than 
the LSDA which includes them roughly in a mean-field 
manner (and in any case the ASA then is a serious issueL 

Recent photoemission investigations of Tsujioka et acl 
found large splittings of 5 eV between the main peaks 
in the occupied and empty Ci-d densities of states. The 



authors compared that number with the corresponding 
splitting in Ref. @ (« 1.5 eV) and in Ref. @ (» 4.5 eV), 
and interpreted that as an indication of strong, impor- 
tant Coulomb correlations. However, this large difference 
is not entirely due to the Coulomb correlations. An im- 
mediate observation that can be .made from comparison 
of the results of Korotin it et alEl and SchwarzEl is that 
the Ref. [7] has not only the unoccupied d-band shifted up, 
but also the weight inside the band redistributed. The 
latter effect can hardly result from the Coulomb U effect, 
because it is present also in our LAPW LSDA density of 
states (Fig. |[ note that splitting here is about 2.5 eV). 

In order to understand the remaining difference, and 
also to explore the sensitivity of the exchange splitting 
to the technical details, we performed a series of LMTO- 
ASA calculations, using first the automatic maximum 
space filling atomic sphere radii, calculated from the 
Hartree potential distribution by the Stuttgart LMTO 
ASA packageo. Our experience of working with this 
package is that this choice nearly always provides much 
better space filling and much smoother crystal poten- 
tials than simple heuristics. In the case of &O2 this 
yielded R Cr = 2.106, R = 1.872, R E1 = 1.728, and 
Re2 = 1-124 a.u. Unlikftr-the setup used in Ref. ^ and 
other ASA calculationaMEll, there were 4 empty spheres 
of the 1st kind and 2 empty spheres of the 2nd kind per 
Cr02- The setup in Ref. p] generated too large sphere 
overlaps, which, [according to the modern understanding 
of the LMTOEjO are above the maximum limit recom- 
mended for a second generation LMTO. However we ran 
the program with that setup as well for the purpose of 
comparison. The results are shown in Table |. We find 
that this setup consistently overestimates the exchange 
splitting by 0.2 eV. Another wl eV then comes from U, 
which may be considered as an estimate of correlation 
effects beyond the LSDA in this system. It is worth not- 
ing that the Generalized Gradient Approximation, dis- 
cussed in more detail later in paper, takes correlations 
into account in a more sophisticated, but still mean-field 
manner, and leads to an exchange splitting about 0.3 eV 
larger than in LSDA, closer to LSDA+U. The gap in the 
spin-down channel, however, is the same in the GGA as 
in the LSDA, unlike LSDA+U. 

Let us now look at the electronic structure of the spin- 
up subband at the Fermi level. Fig. |^ shows the calcu- 
lated Fermi surface. The three bands form a number of 
pockets; the main motif is a chain of alternating hole and 
electron pockets extending along k z direction. The elec- 
tron pocket ( "pseudocube" ) is centered around T, and it 
touches a hole pocket at the point ko (0, 0, OAtt/c). 
These hole pockets extend all the way up to the Z point 
where they touch the same pocket in the next Brillouin 
zone. There are 0.12 electrons per formula unit, and the 
compensating amount of holes. Interestingly, while the 
T-pocket is quasi-isotropic, this hole pocket looks more 
like a pseudosphere with "tentacles'}— directed towards 
the center of the TXMXZRAR prismEl, that is towards 
the point (w/2a, 7r/2a, tv/2c). The tentacle has a com- 
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plex shape and extends from the above mentioned point 
along (110) direction hugging the pseudocube. The whole 
pocket resembles a hammerhead shark with one belly, 
four heads, and no tails. The main contribution to the 
DOS comes from the "heads" . Two additional small con- 
centric electron pockets occur between A and M. These 
two Fermi surfaces are degenerate on the XMAR face 
of the Brillouin zone. The reason for that is the exis- 
tence of a twofold screw axis in the rutile structure; as 
was shown by Herringcj, due to the time reversal sym- 
metry, this leads to degeneracy of electronic bands on 
the plane perpendicular to the screw axis. This degen- 
eracy, generally speaking, is lifted linearly in k, and the 
momentum matrix elements between the bands are fi- 
nite. This leads, in principle, to a finite static interband 
conductivity! 2 ^, cri n t er (ui — > 0) 7^ 0. Since the length 
of the line where these two pockets touch is small, this 
contribution is small, but still calls for additional cau- 
tion when analyzing experimental infrared conductivity 
in terms of a Drude model. 

The calculated DOS, as shown on Fig.|], has a mini- 
mum ("pseudogap"), but the Fermi level is not exactly at 
the minimum; the calculated N(0) = 1.9 states/eV cell 
corresponds to a linear specific heat coefficient 7(0) = 
2.24 mJ/K 2 mole. This value is sensitive to the zone 
sampling andj-care was taken to ensure its convergence. 
Tsujioka et aB have measured specific heat at low tem- 
perature and fitted their data to the following expression: 

7(T) = 2.5 mJ/K 2 + 0.0169mJ/K 4 T 2 + 1.36mJ/K 5/2 Vr, 
the last term being a magnon contribution. It is worth 
noting that this curve, due to the square root term, meets 
the origin with infinite slope, and already at T = 0.1 K, 
for the parameters given in Ref. [| j(T) w 3 mJ/K 2 . 
This at least in part may explain the disagreement be- 
tween various authors about this value. The calculated 
7(0) is in close agreement with the experiment, which is 
surprising (one expects substantial renormalization from 
electron-magnon interactions), but less surprising than 
the anomalous renormalization calculated by Lewis et 
au. The number given by Tsujioka et au may be an un- 
derestimate, given the sharp change of j(T) near T = 0. 
We, however, think that this number is more realistic 
than the 2-3 times larger values quoted in Refs. [15 26 
On the other hand, the above mentioned minimum in 
N(E) occurs at only 30 meV above the Fermi level. The 
strong variation in this region and high sensitivity of the 
exchange splitting to the computational parameters to- 
gether lead to an extraordinary sensitivity of the calcu- 
lated N(Q) to the exchange-correlation potential (partic- 
ularly, to gradient correction). In fact, in GGA calcula- 
tions the Fermi level falls exactly at the minimum and 
the calculated DOS is only 0.95 states/eV cell, half of 
the LSDA value. Depending on the exchange-correlation 
functional used, the calculated bare electron specific heat 
coefficient can be anything from 1.1 to 2.3 mJ/K 2 mole, 
and correspondingly the needed mass renormalization 
anything from 2.5 to 1.1 . 



The pseudogap in the DOS also leads to a minimum in 
the plasma frequency as a function of the Fermi energy 
position. At the actual position of the Fermi-level, the 
plasma frequency is calculated to be 2.1 eVEO. This is 
smaller than the value (3.35 eV) derived by Basov et a£J 
by integrating the measured optical conductivity up to 
12 000 cm -1 (1.49 eV). The reason is that, as mentioned 
above, there are interband transitions with vanishing en- 
ergy in this system, and a large part of the spectral weight 
integrated up to 1.5 eV comes from these interband tran- 
sitions. Note that (Fig.^J) the plasma frequency is much 
less sensitive to the position of the Fermi level than the 
DOS — changes in the value of N (E) are compensated by 
the opposite changes in the Fermi velocity. In the above 
mentioned GGA calculations, despite a drastic change in 
N(0), ftp was only 5 and 7% smaller for the xx and zz 
components, respectively, than in LSDA. 

Combined with the small residual resistivity (20 fifl 
cm in static measurementsEj and 10 fifl cm tpdess when 
extrapolating the optical conductivity to zeralil), the cal- 
culated Lu p calls for an unusually long relaxation time at 
low temperature: 7 = %/t = 50 to 100 cm -1 (the first 
number is based on DC resistivity and the second on the 
optical data). Basov et at3 fitted their low- frequency 
data (u < 150 cm -1 ) by a Drude formula with oj p w 1.7 
eV and 7 « 30 cm -1 . So the LSDA calculations again 
fall within the range of the the different experiments. 

Another characteristic that can be extracted from the 
calculations is the average Fermi velocity. Note that one 
cannot reliably extract this number from the optical ex- 
periments, since the required effective mass is unknown. 
Although the ratio of squared plasma frequency and the 
DOS, formally speaking, gives the average squared Fermi 
velocity, this procedure becomes misleading when the ef- 
fective mass is rapidly changing over the Fermi surface: it 
overweights heavy-mass parts of the FS and reduces the 
calculated (v 2 F ), while the actual transport is dominated 
by the lighter carriers with higher velocity. However, 
even for the lightest band at the Fermi level the calcu- 
lations give ~ 0.2xl0 8 cm/s. The estimation of Ref. 

(0.6-0.9xl0 8 cm/s) is substantially higher; we ascribe 
this partially to the neglect of the effect of interband 
transitions and the deviation of the carrier mass from 
the free electron mass (and indeed 10 8 cm/s is an unusu- 
ally large number for a d-band metal). Correspondingly, 
their estimate for the mean free path at T = 10 K should 
be reduced from 1500 to 300 A, which is not extraor- 
dinary large. The former number is hard to explain by 
freezing out of spin-flip impurity scattering (due to the 
half- metallic character of Cr02), as it is hard to imagine 
that all charge defects are absent and the only possible 
scattering channels are associated with the spin flip. 

On the contrary, the temperature dependence of the 
resistivity indicates a significant role played by spin scat- 
tering at finite temperature. As shown in Ref. [l^, the 
resistivity below the room temperature can be fitted by a 
quadratic dependence, or slightly better with a power law 
with an exponent n = 2.3 — 2.5. Various models for the 
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scattering by spin fluctuations give exponents between 2 
and 3.5, but in any event much larger that one. Interest- 
ingly, p{T) for CrC>2 looks very similar to that of other 
so-called "bad metals". A typical example is SrRuC>3, 
with p(T) closely reminding p(T) in CrC>2, including the 
change of the slope at the Curie temperatureEI. Inter- 
estingly, in Ref . ^ the observed resistivity was argued to 
be not related to magnon scattering because the prefac- 
tor of T 2 seemed too large (~ 10 — 8 cm/K 2 ) compared 
to the corresponding term for elemental fenramagnets 
(~ 10~ n f2 cm/K 2 ). As we discussed elsewhereca, the ori- 
gin is coupling of magnetism and covalent bonding via O 
Stoner factor. The physics of CrC>2, despite its being also 
a bad metal, is very different. What these two materials 
(and probably other bad metals as well) have in common 
is that in both metals changing magnetic order has strong 
qualitative effects on the band structure, and because of 
that electrons are scattered anomalously strongly by spin 
fluctuations (in this case electron-magnon interaction is 
due to what can be called "magneto-covalent" coupling) . 
Note that the T 2 prefactor in Cr02 is of the same order 
as in SrRuC-3 (about three times smaller). 

Optical spectra of Cr02 reveal a number of interest- 
ing featuresLH. Two features observed in the high-energy 
part of the spectrum are a peak at « 3 eV, easily identi- 
fiable with transitions across the spin-minority gap, first 
pointed out by Uspenski et alE3, and a weaker peak (es- 
sentially, a shoulder) at w 2 eV. This latter feature coin- 
cides both in position and relative strength with a peak 
calculated in Ref. ||l] in the polarization E|| c, coming 
from the spin-majorfUL.band (this feature is weaker in 
the other polarizationo). The third observed feature is 
a broad maximum between 500 and 10000 cm -1 (0.06 to 



1.2 eV), interpreted in Ref. [Eyas intraband conductivity 
shaped in a non-Drude way by a frequency-dependent 
relaxation time. Howevex- a closer examination of the 
results of Uspenski et alE3 reveal an interband feature 
starting at about 0.1 eV and continuing to 0.9-1 eV, at 
least in the E_L c polarization. The relative strength 
of these three features in the calculations may be ques- 
tioned, particularly the excessive weight in the lowest 
one. We should keep in mind that a substantially greater 
error is expected in the calculation of the matrix elements 
within LMTO-ASA technique, than in the calculations of 
the band structure itselfiy. Besides, there is the substan- 
tial polarization dependence revealed by the calculations. 
This introduces additional uncertainty in the Kramers- 
Kronig procedural 

Uspenski et alx3 did not discuss the low frequency 
(oj < 500 cm" 1 ) reg-ien in detail. However, following the 
work of Basov et aO, one may look at this region more 
closely. This requires knowledge of the band structure 
near the Fermi level with a better accuracy than that 
provided by the LMTO-ASA method. Besides, LMTO- 
ASA introduces uncontrollable errors into optical matrix 
elements, hence the need for full potential optical cal- 
culations. We have performed such calculations using 
the Wien LAPW_ code with the optical extension, de- 
scribed elsewhereEl The band structure calculated with 
this cade differs from the one calculated using the LAPW 
codeO used in the transport calculations by less than 1 
mRy. 

As mentioned, Cr02 is special in the sense that it has 
non-zero interband conductivity at u> — > 0. To demon- 
strate this, we show in Fig^ the calculated joint density 
of states (JDOS), defined as 



•AjH = ^2 /(£k l )/(-£k J )<5(£ki - £ki ~ M) C 1 ) 

k 

which is proportional to the optical conductivity due to the i —> j interband transitions in the constant matrix element 
approximation: 

2 

5ijM = -7? J2p 2 i3 f(^i)f(-euj)S(s kj - e k i - huj) = —^-J 13 {uj). (2) 



Here V is the volume of the unit cell, and pij is the 
dipole matrix element for given polarization. As one can 
see from Fig.^, there is indeed a linear in us component. 
A closer look uncovers that there is also a quadratic in 
to component in JDOS. This comes from the transitions 
around the point (0,0,0.47r/c) where two bands touch (the 
"pseudocube" and the belly of the "shark"). Even more 
interesting is the steep increase of the JDOS above 700 
cm" 1 . This is the threshold for the transitions from the 
"hammers" into the "pseudocube" . The weak dispersion 



in one direction of the band forming the "hammers" pro- 
vides a large phase space for this transition. 

The calculated optical conductivity including matrix 
elements is shown in Fig.||. In the following discussion 
we shall concentrate at the low energy (u> < 4.5 eV) re- 
gion, studied in Ref. This region is blown up on 
Fig.^, where we also show the experimental data from 
Fig. 2 of the Ref. |ll[ We have added a Drude part corre- 
sponding to the Co = 6000 (f2 cm) -1 and the calculated 
plasma frequency of 2.1 eV. First of all, we observe that 
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all essential features visible in the experiment are repro- 
duced: a narrow intraband peak with no deviation from 
Drude shape, interband transitions within the i2 g mani- 
fold, starting from u> = and extending up to 1.5 eV (1.3 
in the experiment), followed by a wide absorption band 
between 2 and 4.8 eV (between 1.7 and 4.5 in the exper- 
iment), which is associated with transitions between the 
bonding and antibonding e g bands. The e g manifold in- 
cludes some dispersionless bands (which are understood 
in the tight binding model described in Ref. 0). These 
in turn result in an absorption peak near 2.5 eV. It is 
tempting to associate this peak with a shoulder at 2 eV, 
found in the experiment, although by doing so we risk 
indulging in overinterpretation. Nonetheless, in this case 
the broad hump in the experimental curve, centered at 3 
eV, may be associated with the feature in the calculated 
a zz at 3.5-4 eV. It is notable that the experimental fea- 
tures systematically occur at frequencies 10-20% lower 
than the corresponding calculated ones. This must be 
due to many body renormalization effects. Contrary to 
the assumption of Refs. [?] and [| these renormalizations 
seem to be related with the exchange of virtual magnons 
rather than with the Hubbard repulsion, as this would 
force apart occupied and unoccupied states and would 
shift optical transitions to higher energies, (cf. Eliash- 
berg renormalization due to electron-phonon interaction; 
however, because of the high frequency of the spin exci- 
tations, this effect is seen even far away from the Fermi 
level). Again, one should exercise caution here, because 
the crystalline anisotropy, barely noticeable in the intra- 
band conductivity, increases with frequencies (as one can 
see from Fig.^J) . First of all, we do not know with which 
weights the two polarizations appear in the experimental 
curve. Second of all, even if we knew that, and would 
have averaged the two curves in Fig j^, this procedure 
would not have produced the same results as averaging 



and obviously is zero below the spin-flip transition 
threshold. However this interpretation relies on a thresh- 
old of 500, or at most 1000 cm -1 , while our LSDA cal- 
culations give at least 0.25 eV, that is 2000 cm -1 . This 
number is defined by the exchange splitting and the many 
body correlation effects neglected in the LSDA are likely 
to increase it even further (cf. Ref. fj]). Interpretation 
in terms of interband transitions fits the experimental 
picture much better. 

The level of agreement between the observed low- 
frequency plasma energy, 1.7 eV, and the calculated value 
of 2.1 eV is satisfactory. It is tempting to ascribe over- 
estimation of the plasma frequency to many body effects 



two reflectivities and restoring the conductivity from the 
average reflectivity (which would correspond to the ex- 
perimental procedure). To illustrate that, we compare 
the average calculated reflectivity (assuming 1:1 weight- 
ing of R x and R z , that is, an orientation preferential for 
z polarization) with the reflectivity calculated from the 
average dielectric function (Fig. [t]). While in general the 
curves are similar, the latter overemphasizes the feature 
at 1.3 eV, right below the screened plasma edge at 1.6 eV. 
Indeed, we observe from Fig.^j that the structure in the 
experimental curve in this energy region is substantially 
weaker than in the calculated one. 

The last note regarding comparison with the exper- 
iments regards the low-temperature conductivity dis- 
played in Fig. 4 of the Ref. [Ly. As we have noted be- 
fore, the Fermiology of Cr02 is unusually sensitive to the 
approximations used for the exchange-correlation func- 
tional. In particularly, the above mentioned smaller den- 
sity of states in GGA is due to the fact that the ham- 
merheads of the hole pockets entirely disappear when the 
Fermi level is raised by 30 meV. So, one should take the 
calculation in the far infrared (uj < 500 cm -1 w 60 meV) 
region with a grain of salt. With this in mind, the agree- 
ment between the frequency of the threshold (700 cm -1 ) 
seen in FigJJ and the frequency £3 w 300 cm -1 , where 
the experimentally determined low temperature optical 
conductivity starts to deviate from the Drude law, is 
very reasonable. Note that the original interpretation 
of Basov et al was in terms of mfraband transitions and 
a strongly frequency-dependent relaxation time. They 
assumed a threshold in the relaxation time around 500 
cm -1 and ascribed this threshold to the minimum energy 
for spin-flip scattering. If the main scattering mechanism 
is by magnetic impurities, then indeed it is easy to show 
that in the constant scattering amplitude approximation 
the relaxation time depends on frequency as 



(3) 



(cf. high-T c cuprates where the LSDA underestimates 
uj p by a factor of two). However, one should be very 
cautious about such an analogy, because renormalization 
in cuprates affects the quasiparticle velocity but not the 
area of the Fermi surface, so that the DOS (oc J <1Sf/vf) 
is underestimated by about as much as (oc J vfcISf) 
is overestimated. Here, however, both N(Q) and uip seem 
to be overestimated, suggesting that it is the total area of 
the Fermi surface that is overestimated. Unlike the high 
T c cuprates, the Luttinger theorem does not fix the size of 
the Fermi surface because both electron and hole sections 
are present in Cr02- The most natural explanation of the 
overestimated plasma frequency in the LSDA is thus not 



/ui/2 
[Ni(s - w/2)iV T (£ + oj/2) + N^e + w/2)7V T (e - w/2)] ds, 
-U>/2 
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a failure of the one-electron picture, but inaccuracy of 
the density functional approximation used in the calcu- 
lations. The fact that GGA calculations yield smaller 
(substantially) N(Q) and (somewhat) ui p illustrates this 
point. We should emphasize that we do not claim the 
GGA is necessarily a better approximation than LSDA 
for Cr02; it does give larger separation between the oc- 
cupied and unoccupied d states, both in the spin-flip and 
in the same spin channels, thus being closer to LDA+U 
than LDA is. We are, however, sceptical about superi- 
ority of LDA+U in case of Cr02, and indeed the optical 
data suggest that the LSDA bands should be corrected 
in the direction opposite to LDA+U or GGA, that is, the 
band widths should be smaller, as well as the separation 
between occupied and unoccupied bands. 

Finally, an interesting question arises about the evo- 
lution of the optical and transport properties above the 
Curie temperature. One may expect that at sufficiently 
high temperatures the optical response would be deter- 
mined by the paramagnetic bands structure rather than 
by the ferromagnetic one, but this temperature may be 
much higher than Curie temperature. In fact, local mo- 
ment systems, e.g. Fe, retain magnetic splitting at the 
atomic level well above the Curie temperature where the 
spins lose their long range order. One way to detect such 
a situation in zero temperature calculations is to com- 
pare ferro- and antiferromagnetic energies and moment. 
If the difference between the two (set by the interatomic 
exchange energy) is small compared to the intraatomic 
exchange (Hund energy) taken as the difference between 
a magnetic and the nonmagnetic state, and the moments 
are not sensitive to the local order, one has a local mo- 
ment system. 

Cr02 is an intermediate case: the antiferromagnetic 
structure appears higher in energy than the ferromag- 
netic by about a quarter of an eV/cell (which corre- 
sponds to J ~ 200 K, in accord with the experiment), 
while the one without spin polarization is higher again 
in energy thais-the antiferromagnetic one by about the 
same amountEj. On the other hand, the antiferromag- 
netic Cr02 is metallic (not a half- metal), but the mag- 
netic moment inside the Cr atomic sphere is only weakly 
changed (1.56 hb vs. 1.97 (J-b)- 

To summarize, we present the first all-electron and full 
potential LSDA band structure and optical coefficients 
calculations for Cr02- We show that apart from anoma- 
lously strong electron-magnon scattering, which natu- 
rally follows from its half-metallic band structure, the 
experimental observations may be interpretable in terms 
of quasi-one-electron band theory. Thus there is no ex- 
perimental smoking gun in regard to strong correlations 
related exotic phenomena in Cr02- Contrary to some 
conclusions deduced from less accurate band structure 
calculations, we find very modest specific heat renormal- 
ization, and the optical plasma frequency in reasonable 
agreement with the low temperature far infrared opti- 
cal experiments, and interband transitions starting from 
zero frequency. The latter explains the non-Drude shape 



of the observed optical conductivity in the range of 500 
to 4000 cm- 1 . 
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FIG. 1. Calculated band structure of ferromagnetic Cr02 
in the spin-majority channel. The Fermi energy is denoted by 
the dashed horizontal line at 0. 
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FIG. 2. (a) LAPW DOS of ferromagnetic Cr0 2 on a per for- 
mula unit basis. The shaded area is the Cr d contribution. The 
Fermi energy is at 0. (b) Blowup of the DOS for spin-majority 
channel inside the spin-minority gap. Also shown are plasma 
frequencies for the two crystallographic directions as function 
of the chemical potential. 



7 




FIG. 3. LAPW Fermi surface of Cr0 2 . Note the electron 
"pseudocube" around the V point and hole "hammerheads" 
touching the "pseudocube" between T and Z. Small pockets 
along the AM line. 
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FIG. 5. Calculated optical conductivity of Cr02- Spin- up 
and spin-down channels are shown separately. 
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FIG. 6. Calculated optical conductivity of Cr02- Spin- up 
and spin-down channels are added together. Symbols show ex- 
perimental data at 300 K from Ref. [ll| Fig. 2. 




FIG. 7. Polarization average of the calculated resistivity, 
compared with the resistivity corresponding to the polarization 
averaged calculated dielectric function. Averaging over polar- 
ization is performed assuming equal weights of the in-plane and 
out-of-plane polarizations. 
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TABLE I. Spin-flip scattering threshold (distance from the highest occupied spin-up state to the lowest unoccupied spin-down 
state) in eV. HL, VWN, vBH, LM, and PW91 denote, respectively, exchange correlation potentials of Hedin-Lundquist, 
Vosko-Wilk-Nussair, von Barth-Hedin, Langreth-Mehl, and Perdew-Wang. The last two are GGA functionals. 



LAPW (HL) 


LAPW (PW91) 


LMTO (vBH) 


LMTO (VWN) 


LMTO (LM) 


LMTO(PW91) 


LMTO (VWN)" 


LMTO(PW91) a 


0.31 


0.69 


0.24 


0.31 


0.42 


0.55 


0.50 


0.62 



a Using setup from Ref.l 
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